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Abstract 

In this paper we tried a different approach to work out the integrals of e x and 
e~ x ™ . Integration by parts shows a nice pattern which can be reduced to a form 
of series. We have shown both the indefinite and definite integrals of the functions 
mentioned along with some essential properties e.g. conditions of convergence of 
the series. Further more, we used the integrals in form of series to find out series 
solution of differential equations of the form — (n— 1) ^ — n 2 x 2n ~ 1 y — nx n = 

and — (n— l)^f— n 2 x 2n ~ 1 y+(n— 1) = 0, using some non standard method. We 
introduced modified Normal distribution incorporating some properties derived 
from the above integrals and defined a generalized version of Skewness and Kur- 
tosis. Finally we extended Starling's approximation to lim (2n)! ~ 2nV2n (^) " 

1. Introduction 

Integrals of exponential functions like e x " or e~ x " can generally be expressed in 
summation form V = — S ppr + k\ and (— l) m — ft ttt + ko respectively. We 

^— ' m!(nm+l 1 ^— ' v ' m!(nm+l) z J 

m=0 m=0 

tried a different approach to evaluate the integrals which can be applied to solve some 
special kind of mathematical problems. When we tried to integrate the functions 
mentioned, by the method of integration by parts, we see a nice pattern which can be 
reduced to a form of a series. 
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India for his invaluable suggestions. S.P. wants to thank "Deutsche Forschungsgemeinschaft" (German 
Research Foundation) [DFG](GK 1177) for funding. 
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In this paper we will show. 
• For all positive integer values of n, indefinite integral of e x " and e~ x " is 



£ 



{-l) r n r x 1+rl ' 
~0 f[ (l+pn) 

p = 



+ Ci and e 



r=0 II 

p=0 



+ C2 respectively [ Vx; at 



least Vra G Z+]. 



• We showed the convergence of the above series to work out the proper integrals 
(e.g. 1.2.1) and shown that the modified interval of convergence for improper 
integrals are ] — 00, oof for n=even. And — 00 ^ x < 00 for n=odd for J e x " dx 
and [—00, 00] for n=even and — 00 < x ^ 00 for n=odd for J e~ x " dx. 



We applied these to solve some special kind of differential equation of the form 
^-{n-\)f x -n 2 x 2n - x y-nx n = and ig-(n-l)|-n¥"- 1 y + 
(n — 1) = as series using some non standard method. 

Further we used some properties of the above integrals to define some proba- 
bility distribution function (P.D.F) and finally we defined a generalized version 
of Skewness, Kurtosis and Kurtosis excess for generalized "(f)*' 1 order Normal 
distribution" i.e. y(x) = p e for even n, 



Coefficient of Skewness 



- Kurtosis = ^ 

mi 



Kurtosis excess = Z ?S L — (1 + n). 



"<:, 



Where m s denotes s t h order central moment of the distribution. 
• We also extended Starling's approximation to lim (2n)! ~ 2n\/27r (^) 2n . 

2. Integration of e x " and e~ x ™ 

Generally integrals of e x " and e _x " can be expressed in an easy summation form 
£ + kl and S (-l) m 4^+I) + fc2 -( 2 ) respectively.But we 

m=0 m=0 

tried to evaluate them in a different manner, which is rather difficult but have greater 
implacability in some special kind of mathematical problems. 

2.1. Indefinite integral 

For all positive integer values of n 

f e x "dx 
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e x "x — n J e x " .x n dx 

e xn x - e xn ^x n+1 + +Jj Je xn .x 2n dx 

e x n rr _ „i° n „n+l _|_ p x n n 2 „2n+l n 3 r x" T 3nj„ 

e x e -he ( n+1) ( 2n +i) x " (n+i)(2n+i) J e - x ax 



= e^x - e xn ^x n+1 + 



A :in ■ i 



(n+l)(2n+l)' 

p i" nf T 3n+1 j n* f „i™ „,4nj 

e (n+l)(2n+l)(3n+l) ^ (n+l)(2n+l)(3n+l) J e >X " X 



E 

Lr=0 



(~l) r n r x 1 + nr 

(l+0n)(l+lra)(l+2n)...(l+rn) 



+ Ci 



g (-l) r n r i 1 +"- 
r=0 fl (1+P»i) 



+ Ci 



Vx; at/east Vn G Z+ (3) 



Similarly: 

J e~ xn dx 



X + ^ (l+0n)(l+ln)(l+2n)...(l+rn) 



+ c 2 



E -£ 



=0 11 (l+pn) 

P=0 



+ c 2 



Vx; atleast Vn G Z + (4) 



Note that equation no. (3) and (4) are valid for all x for all n other than[— 1 ^ n ^ 0] 



2.1.1. Observation: 

Taking n= odd and putting x= -y i.e. dx= -dy in equation no. (3) we get, 

Je(-yr d (- y ) 



E 



(-l)-(-l)i+"-nV+" 



=0 n (i+p«) 

p=0 



g (_l)l + (r t + l)r- nrj/ l + nr 

r=0 fl (l+pn) 

p=0 



+ Ci 



+ Ci 



= -p-v 



E 



r=0 II (1+P«) 

p = 



+ Ci 



or, J* e y " dy = e 



r=0 I! (l+pn) 

p = 



Ci 



as n is assumed to be odd 
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or, J e x dx = e 



E 



=0 11 (l+pn) 

p=0 



Ci 



(4A) 



Comparing equation no. (4) and (4A) we have Ci = — C2 



2.2. Definite integral 
2.2.1. Proper integrals: 



(-l) r n'V+ T "' 



E 

r=0 II (1+pn) 

p=0 



For proper integrals, convergence of e 
finite values of x, are essential and convergence of e 



E 



and e^" 



(-l) r n'-a: 1 - 1 



E 



r=0 I! (1+pn) 

p=0 



for 



r=0 I! (1+pn) 

p = 



depends upon 



both e x " and 



j2 (-iyn r x 1+n - 

r=0 EI (1+pn) 

p=0 



We thus performed the ratio test of both e x " and 



E 



=0 11 (l+pn) 

p=0 



which shows 






lim 

r^oo 




= lim 

r— »oo 






x"»-(r+l)! 



lim ^ \x n \ = 0(< 1) 



= lim 
= lim 

r— +00 



ttr+l 



lim 

r— >oo 



p = o 



n r a; 1 +"'" J] (1+pn) 



1+rt+rn 

= lim 1 , w , |-x n | = 0(< 1) 
So, it converges for all finite values of x. 



Hence, e a 



(-l) r n r x 1+rl 



r=0 I! (1+pn) 

p = 



E 



i.e. J e x " dx converges for all finite values of x. 



Similarly it can be shown that f e x "dx also converges for all finite values of x. Same 
conclusion can be drawn from integrals of form of (1) or (2) by ratio tests. 



2.2.2. Modification of interval of convergence: improper integral 

Since the functions e x " and e~ x " are having no point of singularity, finite integrals are 
improper when any (or both) of the limits are infinite. 
Now, 

+00 +00 

J e x "dx= J e x "dx+ J e x "dx (5) 
—00 —00 
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oo 

Now, by n test of J e x " dx....(6) we see lim x^e x " = oo for, \x = or 1(<1) (where 

x assumes only positive value) is evidently diverging, 
o 

Again, for f e x " dx ; x Assumes only negative values in this interval 

— oo 

Let, x=-y; i.e. dx=-dy (y > 0) 
2.2.2.1. for n = even : 

oooo ( b b ) 

J e x "dx = J e yn (—dy) = J e y " dy = J e x "dx < As, J f(y)dy = J f(x)dx> 

— oo oo ! a a J 

By earlier conclusion it diverges 



2.2.2.2. for n = odd 

oo oo 

J e x " dx = J e~ yn (—dy) = j e~ y "dy = J e~ x "dx. 

— oo oo 



By j! test it is evident that 



lim x^e 

x— >oo 



lim 

rr^oo 



For n = 2(> 1) 



lim 



As,jf(y)dy = jf(x)dx 

a a 



(7) 



For n < ,It is oo ( diverges) (n need not be integer even). But we have not 
considered so. For n > 0, it can be shown (not necessary that n is odd, it may be 
even or fraction also) that expression (7 ) is convergent For integer value of n, the 
expression obviously is 0, hence convergent. For fractions V varies from to oo so, 
what ever the value of n may be, as it is finite, there will be r's > ^. So, 



lim 

rr^oo 



lim 

x— >oo 



-2 , 'In 1 



• ^yya; +oo+oo+c^ 



= 



(7A) 



Hence convergent. 

oo 

So, J e~ x "dx ,is convergent Vn G Z + i.e. for odd n's also. 



i.e. f e x "dx converges. 

— oo 

So modified interval of convergence of f e x " dx is ] — oo, oof for n = even. And 



— oo ^ x < oo for n=odd. 

+oo +oo 

Now J e~ x "dx= J e~ x "dx+ J e~ x "dx 

— oo — oo 



(8) 



As shown earlier, J e~ x "dx is convergent Vn G Z + and for, f e~ x " dx, x Assumes 

oo 

only negative values in this interval. 
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Let, x=-y; i.e. dx=-dy (y>0) 



2.2.2.3. for n = even : 

oo oo 

/ e~ x "dx = J e~ y "(—dy) = J e _2/ " dy = J e~ x "dx. 

— oo oo 



b b 

AsJ f(y)dy = f f(x)dx 

a a 



By earlier proof it is convergent. 



2.2.2.4. for n = odd : 

oooo C b b } 

f e~ xn dx = f ey n {-dy) = f dy = f e x "dx < As J f{y)dy = f f{x)dx \ 

— oo oo la a J 

Which diverges by earlier proof. So modified interval of convergence of f e~ x " dx is 
[—00,00] for n= even and — 00 < x ^ 00 for n=odd. 



2.2.3. Observation 



E 



(-l) r n r x (1+nr:i 



=0 n 

p=0 



E x nr 

r=0 



r=0 fl (1+pn) 

p = 



1 4. *1 
1 + 1! 



Hr" + + 



+ 



(l+0n) (l+0n)(l+ln) ^ (l+0n)(l+ln)(l+2n) 



-...OO 



+ 



0!(l+0n) 



„l+0n 



0!(l+0n)(l+ln) l!(l+0n) 

i n 

+ 



x 



l+ln 



0!(l+0n)(l+ln)(l+2n) l!(l+0n)(l+ln) ^ 2!(l+0n) 



X 



l+2n 



.OO 



= E (-1)" 

m=0 

Similarly 



E 



(-1)""" 



r=0 r! (1+pn) 

p = 



,1+mn 



E -£ 



r=0 II 

p=0 



= E 



m=0 



E 



(-l) T, n m - r 



r=0 r! (1+pn) 

p=o 



.T 



1+mn 



From equation no.(l)and(9), comparing coefficients of x 1+mn we get, 



E 



(-l) p n" 



1 = r_iy 

m!(nm+l) \ / 

From equation no. (2)and(10) comparing coefficients of x 1+mn we get, 



r=0 r! n (1+pn) 

p = 



\m 1 



E 



(-i) r w 



C_lf 

V / m!(ran+l) 

Obviously equation (11) gives equation (12) if multiplied by (— l) m . 



r=o r\ n (i+pn) 

p = 



(9) 



(10) 



(11) 



(12) 
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3. Application of integrals of e x " and e x " in the form of series as series 
solution of some differential equations. 

Considering two second order differential equations: 

*§f -(»-!)£- n 2 x 2n ~ x y - nx n = 
or, 

g-^S-n 2 * 2 - 2 ^^- 1 (13) 
and 

x & - ( n ~ l )t - nV"-^ + (n - 1) = o 

or, 

S-^^-^ 2 "^ = ^ (14) 

Here equation no. (13), (14) are actually the special forms of the equation 
+ P^ + Qy = R (15) Where P,Q,R are functions of x only. 

By method of change of independent variable and taking some arbitrary function 
z = J e~~ I Pdx dx equation no. (15) can be reduced to ^| + Q±y = R\ (15-A) 

where, Q 1 = R x = -g-? 

V dx J V dx J 

Now, for both (13) and (14) z = ^ i.e. x = (nz)* = let, / _1 ( z )- 
Now, (§) 2 = x 2n ~ 2 or, Qi = -n 2 

So, (13) becomes ^| — y = R\ (16) 
Now R\ = nx 1 ~ n = rj w"m„-i for equation (13) and = ,, for equation 

(14) 

So, equation (13) becomes 

§-^ = Fl- ^ 

and equation (14) becomes 

#v _ _ -(n-l) /-, H) 

Evidently (17) and (18) are second order linear differential equation with con- 
stant coefficients. The complete solutions of these are given by Complimentary func- 
tion^ .F . ) +Particular integr al(P.I . ) 

For both, C.F. is solution of the equations. 

- n 2 y = or, symbolically, (D 2 - n 2 )y = {D = £} 

or, C.F.= k\e nz + A^e - ™ 2 = kie x " + k2e~ x " ki, &2, arbitrary, constants. 



Now P.I. of equation no. ...(17) is (P.L)i = miz^i) 



and for equation no. ...(18), (P.I.)?. 



l 

(D 2 -n 2 ) 



-("-!) 



{f-Hz)V n - 



It seems to be complicated to be evaluated in this approach. So alternative 
methods can be taken. By definition particular integral is any particular solution of 
the given differential equation. So (P.L)i is a particular solution of (17) or (13), where 
as (P.I.)2 is any particular solution of equation no. (18) or (14). 
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It can be shown that g(x) 
equation no. (13). 



E 



=0,2,4,. .oo EI (1+pn) 

p = 



is a particular solution of 



Where as, f(x) 



E 



r=l,3,5,..oo II (1+pn) 

p = 



is a particular solution of equation (14). 



By equation no. (3) and (4) 

/ e x "dx = e x " [f(x) - g(x)} + d 

and 

Je'^dx = e~ x " [f{x) + g{x)] + C 2 
Differentiating equation (19) with respect to x, 
e x n _ e x n [y'( x ) — g'{x)] + e x " [f(x) — g(x)] nx n ~ x 
Differentiating equation (20) with respect to x, 

e -x» = e -x« + g l( x ^ _ g-as" + 5 ( x )] nx n-l 

or 

[/'(x) - </(*)] + [f(x) - g(x)} nx n ~ x = 1 
and 

\f{x) + g'{x)\ - [f(x) + 5 (x)] nx n ~ l = 1 
from equation (23) & (24) 
/'(x) = nx n ~ 1 g(x) + 1 
</(x) = nx 1l ~ 1 f(x) 

Differentiating equation (25) & (26) and using (25) and (26) we have, 
g"(x) - ^g'ix) - n 2 x 2n - 2 g{x) = nx^ 1 
and 

fl{x) _(n-l) fl{x) _ n 2 x 2n-2 f{x) 



-("-!) 



So, 

S(z) 



E 



r=0,2,4,..oo II 

p=0 



is a particular solution of (13) 



where as, 



/(*) 



E 



r=l,3,5,..oo II (1+pn) 

p=0 



is a particular solution of (14). 



Now the complete solution of (13) is y = k\e x " +k2e x " + 

And that of (14) is y = k\e x " + A^e -2 ^ + 
Where, k\ and k2 are arbitrary constants. 



(19) 

(20) 

(21) 

(22) 

(23) 

(24) 

(25) 
(26) 

(27) 
(28) 



E 



r=0,2,4,..oo il (1+pn) 

p=0 



r=l,3,5,..oo fl (1+pn) 



{28A) 
(28B) 
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3.1. Graphical representation of y(x) = e _I " for even values of n. 




1 — 1 — i — ' — i — 1 — i — ' — i — ' — i — 1 — i 

-6-4 -2 2 4 6 

X 

Figure 1: Values of y(x) for n= 2,4,6 

Fig 1. shows the curves of y(x) = for different even values of n. 
Again, we have, 

When n tends to oo, in that case 

• y(x) = lim e~ x " = 1, for - 1< x < 1 

n— +oo 

• = \, at,x = ±1 

• = 0, elsewhere 

Which describes a rectangular function [fig 2.]. 
But if n is very big(yet finite), 
lim e~ x ™ = lim e~ x " = lim e _a; " = - 

x^±l+ x-»±l- x-t-±l e 

To be noted that taking a big value of even n the aforesaid rectangular function 
can be approximately replicated in the field of computation. The value of n is to be 
chosen depending upon the ability of the computing machine. 

4. Admissibility of y(x) = Ae~^ ° ) as a Probability Density Function 
(P.D.F) in (— oo < x < oo) for positive, even values of n. 

_ J_ / x-m \ n 

Let us, consider a function y(x) = Ae -I ° > (A being constant) y(x) must 
hold the following properties in order to be a Probability Density Function (P.D.F) in 
(a,b). 

oo 

I. y(x)>0 in(a,b), II. y(x)=0 else where and III. f y(x)dx = 1 

— oo 
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2 -, 






>- 

4- 




(-1,1 le) \ 
, \- 


B- 


• (1,1/e) 

i 1 , 



Figure 2: y{x) = lim e 



It is known that normal distribution is a bell shaped continuous probability dis- 
tribution and it must hold the following properties. 

• Normal curve has two parameters, mean = median = mode = 'm' and standard 
deviation a. 



normal curve is of form y(x) 



l 



OV27T 



e ^ 



(29). 



• the curve is symmetrical about x = m. 

• odd order central moments are 0. 

• Coefficient of Skewness and Kurtosis excess are 0. 



Here, the Probability density function in form of standardized variable z = 

is = 7fc e " f2 (3 ^- 

For y(x) = Ae~™( ° ) , (n being even) conditions I and II are automatically 
satisfied. In order to satisfy condition III, constant A is to be chosen suitably. So 
that, 

!- 00 y{x)dx = i 



roc 

-oo 
or, 

OO 



f Ae-ni^) dx = 1 

— oo 

or, 

— oo 

+oo 

For even values of n, it has been shown that the improper integral J e _x " dx 

— oo 

exists and that can be expressed as 2 T (1 + \).... (30 A). Let that be denoted as P n . 
So, 
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A^a J e -te)" d(^M\ = A^iaP n 

— oo 

or, ^ yfnoP n = 1 or, i = y£^p n 
Therefore, 

= V J^ uP e"( g ) is a Probability density function in (—00 < x < oo)(31) 
Which can be called j th , order (where n=2j) Normal Distribution. Equation (29) 
is a I s * order Normal Distribution in this regard. 



4.1. Properties of distribution corresponding to P.D.F y(x) 



fnaP n 



1 ( x-m \ n 
n\ a I 



Equation no. ..(31) can be further reduced in form of standardized variable z 



is y(z) 



( — 00 < z < 00) 



(32). 



4.1.1. The distribution is symmetric about x=m i.e. z=0 and is bell shaped. 



From equation (32) we have, 

y(z) = 



tynPn 



-OO < Z < OO) 



For even values of n, this is an even function. So, this is symmetric about z=0. 
Which implies equation (32) is symmetric about x=m. 
Again, we have, 



-e n 
1 



So, ^ = 0, at [z = 0,and,y = 0] 

at z = 0, y(z) = n/ l p and at y = 0, e~« = or, z = ±00 



(n = even) 



i.e. 



dz 



0, at 
Now, 

d 2 y(z) = dj-z^y) 
dz 2 _ dz 
■a-\dy 



(°'Ww)]' and ' [( ±0 °'°)] 



n-2 



~n— 2 I y n 



f z +y(n-l)z 

~z n - x {-z n - l y) + y{n - l)z n - 2 } 
~(z n -(n-l))y 



(33) 



Obviously, this expression approaches to from negative side of number line. 
Let e be a variable, which tends to z or, 



d 2 y(z) 

dz 2 



lim U 2n ~ 2 ~ (n - 1) e n ~ 2 ] 



Now, n > 2, or, (n — 1) > 1 , 2n-2 > n-2 and e is a fraction which tends to zero, 
lim [e 2 ™" 2 - (n - 1) e n ~ 2 } < 

or [^Ml < n 

L J (°'^) 

or, y is maximum at z=0 and is ^ p . 
Again, 

Let , 5 1 be a variable, which tends to z and 7 is a variable tending to y. 
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or, 



d 2 y{z) 

dz 2 



= lim lim 7 5 n ~ 2 \5 n -in- 1)1 

(±oo,0) <5^±oo7-+0 



Now, n is even number and > 2 and '5' is a number which tends to ±oo. 



or, 



<Py(z) 
dz 2 



= lim lim 7 \5 2n - 2 - (n - 1) <5 n - 2 l 

(±oo,0) 5^±oo7^0 L J 



(34) 



Obviously expression (34) approaches from positive side of number line, as 7 is 
non negative as it represents y (evident from the expression (32)). 



or, 



d 2 y(z) 

dz 2 



> 



(±00,0) 

So, y is minima at z= ±00 and is 0. 

Therefore, at ^0, ^ p ^ the function gets maxima and at (±00, 0) the function 
gets minima. 

Now equation (33) =0 at points of inflexions. The condition is satisfied for z=0, y=0. 
or, [z n — (n — 1)] = 0. At first two conditions, maxima and minima lie. So condition 
for point inflexion is z n - (n - 1) = or, (z 2 ) 2 - (n - 1) = (34A). 

^ is whole number as n is even. So, equation (34A) may have ^ numbers of 
possible solutions, but at least one possible real solution is. 
(z 2 ) = (n — 1) " or,z = ± (n — 1) » 

there, y = e "{ =t ^ n = e - "^ - ^ = 

So, the curve has two points of inflexions, one in first quadrant {{n — 1)«, \e^j 
another is in second quadrant ^— (n — ^e«^. (Where,n is even no). 



1.5- 



0.0- 



• (n-D ,,n 


(abcissa) 


- (1/e)e 1,n 


(ordinate) 


• -(n-D ,,n 


(abcissa) 



3S -1.0- • 



20 



40 



60 



100 



Figure 1: Abscissas and ordinates for different values of n of 'points of inflexion'. 
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Fig. 3 shows variations of abscissas and ordinates of 'points of inflexion' for differ- 
ent values of n. 
Again, 

dy ^ = —z n ~ l y, n-1 is an odd no; z varies from — oo to oo and y is non- 
negative, [evident from expression (32)]. So the gradient of the curve of P.D.F [ex- 
pression (32)] is positive when z is negative and is negative when z is positive. 

With the aforesaid arguments it is clear that the distribution is symmetric about 
x=m i.e. z=0 and is bell shaped. 

4.1.2. Moments about mean. 

• Odd order moments about mean. 

m 2p+ i = J {x-m) 2 P +l e-^-^) dx 

— oo 

oo 

= _1 f a 2 P +l z 2p+l e -r a dz 
p n una J 
— oo 

Integrand is an odd function. So, all odd order moments about mean are 0. 

• Even order moments about mean, 

We worked out analytically the expressions of even order moments as follows.. 

OO l(x-m\ n 

m 2p = ^h^ J (x-m) 2 Pe-^ — > dx 

— CO 

oo „ 

= — ^ f a 2 Pz 2 P e ~^adz = z\ 

— oo 
oo 



= / z^e'^dz (35 A) 

Pn n — oo 

As both 'n' and '2p' are even so integrand becomes even function. So (35A) can 
alternatively written as.. 

-^Jz 2 Pe-^dz (355) 



00 



The integral J z 2p e « dz exists for any value of 'p'. 







Because, (by fi test) lim z^z 2p e « tends to zero for / u=2(>l).So integral as 

z^oo 

stated in (35A) or (35B) exists for any value of 'p'. 
Now (35A) can be written as.. 



_ z 2p-n-l e - 



—00 
00 



+ (2p-n-l) / z 2p - n e-^dz 



Pn yfn 

= ^^(2p-^T) / a 2 P~ n z 2 P- n e-^dz 



X A11 the graphs are plotted using Origin 7 software package. 
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= (cr n ) (2p+l- n) m 2p - n 
By repeated reduction.. 

m 2p = (cr™) 2 (2p + 1 - n) (2p + 1 - 2n) m 2p _ 2n 

= (cr™) 3 (2p + 1 - n) (2p + 1 - 2n) (2p + 1 - 3n) m 2p _ 3n 

= (cr™) 4 (2p + 1 - n) (2p + 1 - 2n) (2p + 1 - 3n) (2p + 1 - An) m 2p _ in ..etc (36) 
For some positive integer value of 'k' if '2p-kn' assumes zero value, then only by 
repeated reduction the expression (36) gets reduced to 

2p/n 

m 2p = ( a n ) 2p/n n (2p + 1 - kn) m (36 A) 

k=i 

Obviously, [mo = 1] 

If there be no such integer 'k' for which '2p-kn' becomes 'zero', then following 
things happen. Let 'q' be the greatest integer for which '2p-qn > 0'. '2p-qn' must be 
lesser than 'n', else another 'n' could be subtracted. As 'p, q, n' all are integers and 
'n' is even so '2p-qn' is an even integer. 

In this case the ultimately reduced expression 

m 2p = (cr™) 9 f[ (2p + 1 - kn) m 2p . qn (365) 

k=l 

For a given 'n' and 'p' , 2p-qn can assume any even value at most up to n-2. So 
expression (36B) involves terms mo(for qn=2p), m 2 ,m^,mQ m n - 2 

Obviously these can't be further reduced. For a given 'n' , any even moment can 
be expressed in term of one of these as an ultimately reduced form. So these may 
be called fundamental moments". For a given even number 'n', there should be 'n/2' 
numbers of such fundamental even order moments. In that respect 'mo'is the first 
fundamental moment" with value equal to unity. If for a given 'n' and '2p' , 'mo'is 
achieved by successive reduction then by (36A), 

2 P /n 

m 2p = (a n ) 2p/n n (2p + 1 - kn) (36C) 

To be noted that, as the order of target even order moment is increased, any 
particular fundamental moment appears periodically with a period of 'n'. For a given 
'n', 

mo is achievable in [m n , m 2n , m^ n , m^ n ...\ etc. 

m 2 is achievable in [m 2+n , m 2+2n , m 2+ 3 n , m 2+ 4 n ....] etc. 

m4 is achievable in [m4 +n , m4 + 2 ni ™4+3n, m4+4n--] etc. 

similarly, 

[m n _ 2 ] is achievable in [m i ^ +n , m^ +2n , m^ +3n , m^ +An ...] 

If for different values of 'n', m„, m 2ra , ms n , m± n .... etc (i.e. moments involving 
"first fundamental moment" ) are evaluated by reduction relation, then they appear in 
a regular pattern as shown below.. 

n = 2, m 2 = (cr 2 ) 1 -!; . ™4 = (a 2 ) 2 . 1.3; , m 6 = (a 2 ) 3 . 1.3. 5. etc. 

n = 4, m 4 = Or 4 ) 1 -!; , m s = (o" 4 ) 2 .1.5; , m 12 = (o- 4 ) 3 .1.5. 9. etc. 

n = 6, m 6 = (a 6 y.l; , m 12 = (cr 6 ) 2 . 1.7; , m 18 = (cr 2 ) 3 . 1.7.13. etc 

Thus for a given 'n', 
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m n = {a 71 ) 1 . 1- , m 2n = (a™) 2 . 1.(1 + n); , m 3n = (a n ) 3 .l.(l + n).(l + 2n) Thus, 
m fcn = {a n ) k n (1 + rn) (36D) 

r=0 

Even order moments can also be expressed in following forms.. 
By (35B) 




Now from (30A) Pn = %T (£) 



"T-2p = ( V 7 "- ") 



2 P r(i- 



(38) 



4.2. Skewness and Kurtosis. 

Skewness is the measure of degree of asymmetry and kurtosis is the measure of degree 
of peakedness of any distribution, i.e. these two are two parameters telling about the 
geometry of a distribution. For Normal distribution we know that Coefficient of 

Skewness = \ ■•■•[1] 

kurtosis = ^ and kurtosis excess = ^ — 3 [31 

However Pearson expressed Kurtosis as — 3. ...[2] 

Now, Normal distribution is special form of y(x) = - — m) 2p e 

when 'n=2'. For that Skewness is zero as third order (odd order) central moment is 
zero. In that case kurtosis is '3' and kurtosis excess is 0. By (38) kurtosis could be 
evaluates as, 
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(ifna) w^wi 



n ) \ /; / " \ II , 



r(5)r(i) r(3 + 3) /?(i,5) 

- r (3) r (3V r (5 + n " 0/3 3) 
Vra/ Vn/ \n nJ ' Vn'n/ 



(39) 



(40) 



And kurtosis excess= ^4fj^4|j — 3 (4CL4) 

According to this the distribution is always mesokurtic, except for 'n=2'. 

To be noted that for Normal distribution (n=2), both 777,4 and 771,2 are reducible 
in terms of 'first fundamental moment' i.e. 'mo'(=l). But for 'n=4', 'mo'and for 'n > 
4' neither is expressible in term of first fundamental moment. Hitherto it is assumed 
that Normal distribution is the basis of measurement of Skewness, kurtosis or kurtosis 
excess of all distributions of practical use. But keeping the fact in view that, each 
member of the 'family of curves' represented by y(x,n) = - ^_^ e~"( ° ) , ('n' 
is even number and is the parameter of the family) shows bell shaped distribution 
as that of the Normal distribution (the member for n=2). In that respect, to keep 
the 'analogy', it is rational to express kurtosis, Kurtosis excess and Coefficient of 
Skewness in terms of two consecutive even order moments involving 'first fundamental 
moment' i.e. 'too'(=1) for any given 'n'. In that case Kurtosis, Kurtosis excess or 
Skewness values of any distribution obtained from any experiment subjected to fit 
normal distribution will depend upon the order of the normal distribution chosen. 
Then kurtosis of normal distribution of order ' | ' is r ^ L . 

(a") 2 n (1+rn) 

By (36D) ^ = ^ v = 1 + n (41) 

|K) n o a+"0j 

Or, 

by (38) 



m 2n / _/— \2n T + ^) J , , n T + ^) 



2 



m- 



r(2 + ^)-r(i) _(i + i).i 



r(i + ^)-r(i + i) 



1 1 

n ' n 




In that respect Kurtosis excess of any obtained distribution is — (1 + n),and 
that of the chosen normal distribution itself is '0'. 

Coefficient of Skewness of any distribution of form y(x,n) = — l^e"( ° ) 
is 'zero', as the curves are symmetric about the mean of the distribution. Although 
all odd order moments are zero for y(x,n) = — ^^e _ «( ° ) it will be rational to 

v ' Pn \yna 

express Coefficient of Skewness for any given value of n as, , in order to keep 

the 'analogy' and dimensional similarity (of same order). 



INTEGRATION OF E x AND E~ x IN FORMS OF SERIES 



17 



4.3. Multivariate distribution: 



Consider a multivariate distribution for variates z±, Z2, z%, ...z r . with no correlation, 
where Zi = Xi ~ mi . the corresponding orders being m, ri2, 113, ...n r .then combined 
P.D.F is given by 



,-s($) 



n n 2 n 3 



(42) 



For r =2 it becomes a bi-variate distribution Yo 



,"1 "2 



n\ = ri2 = 2, it becomes bi-variate Normal distribution, i.e., Y2 = o^e 2( z ?+ z i)(43) 



5. Extension of Sterling's approximation. 

John Wallis in 1656 [4] has given an expression as follows, 



2.4.6.8.10.. ..00 
1.3.5. 7.9.. ..00 

(2n)M 



n 2 



7T 

2 



or, lim 

n— »oo 



or, lim 



(2n- 1)!! 
(2n)!! 



7T 

2 



n^oo L(2n - !)!!_ 



(44) 



or, lim (2n)!! = lim «/ — (2n — 1)!! 

n— +00 y 2 

or, lim [(2n)!!] 2 = lim A /J(2n - l)!!(2n)!! (44A) 

n^oo n^oo y 2 



Again we have, 



(2n)! = 1.2.3.4.5.6.8 2n 

= (1.3.5.7.9.. ..2n- 1) (2.4.6.8.10.. „2n) 
= (2n-l)!!(2n)!! 

or,(2n- l)!!(2n)!! = (2n)! (45) 
and, 

(2n)M = 2.4.6.8.10.. ..2n 
= 2 n (1.2.3.4.5....n) 

= 2 n n! (46) 
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By the equations (44A), (45), (46)... 

lim \2 n nf = lim J-(2n)\ 

or, lim 2 2n (n!) 2 = lim J-(2n)\ 

n— +oo n^oo y 2 



or, lim (2n)! = lim 2 

n— >oo ra— >oo 



,2n 



or, lim (2n)! = lim 2 

n^oo n^oo 



2n 



(n!) 2 

( lim n!^ 



(47) 



From Sterling's approximation [5] 
By (47) and (48) 



lim n\ ~ y/2im (—) 



(48) 



lim (2n)! ~ 2 2n 

n— >oo 

or, lim (2n)! ~ 2 2 ™ 
or, lim (2n)! ~ 



27m 



(=)■ 



/2n\ 2n 
2m i ( — j 



- -,'2n\ 2n 
or, lim (2n)! ~ 2nv2-7r ( — I 



(49) 



Again in Sterling's approximation (48), instead of 'n' if we put '2n', then, 

f2n\ 2n 
lim (2n)! ~ \/2vr.2n — 

2n-KX> V e / 



. . /2n \ 2n 

or, lim (2n)! ~ V2nV27r — 
2n-KX> V e / 



(50) 



(49)and (50) are two different expressions. Thus Sterlings approximation can be ex- 
tended to (49) in order to approximate the value of lim (2n)\. 
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